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Abstract 

It is a well-known fact that the first and last non-trivial coefficients 
of the characteristic polynomial of a linear operator are respectively its 
trace and its determinant. This work shows how to compute recursively 
all the coefficients as polynomial functions in the traces of the successive 
powers of the operator. With the aid of Cayley-Hamilton's theorem the 
trace formulas provide a rational formula for the resolvent kernel and an 
operator-valued null identity for each finite dimension of the underlying 
vector space. The 4-dimensional resolvent formula allows an algebraic 
solution of the inverse metric problem in general relativity. 

PACS: 02.10, 02.30.T, 04.20. 

I. Introduction 

The eigenvalue problem arises in a variety of different branches of math- 
ematical physics. For instance, it is well-known that quantum systems reach 
stationary states that are the eigenvectors of a suitable linear operator denned 
throughout a complex vector space representing the physical quantum states. 

In the infinite dimensional case real and complex analytical methods have 
been developed to compute the eigenvalues of a linear operator. We mention, 
for instance, global variational methods, which consist of investigating station- 
ary levels of a suitable energy functional, and local perturbative methods by 
means of complex analytic continuation. 

On the other hand, if the underlying vector space is finite dimensional 
the eigenvalue problem becomes an algebraically well-posed problem by means 
of the characteristic polynomial. It is exactly this feature that justifies the 
algebraic approach in this work. Following this general device, the algebraic 
environment is set by the complex number field. After working out a few 
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tools, the most important of these being Newton's identities, and setting very 
fundamental concepts and basic results in complex linear algebra, we succeed in 
achieving the mathematical core of this work: a recursive algorithm to compute 
the coefficients of the characteristic polynomial as algebraic functions in the 
traces of the successive powers of the linear operator. 

With the aid of Cayley-Hamilton's theorem the trace formulas already ob- 
tained provide an operator-valued null identity for each finite dimension of the 
underlying vector space. As a by-product it sheds light on the algebraic struc- 
ture of the associative algebra of complex linear operators. The computational 
skill just attained is enough to yield an operator-valued polynomial with ra- 
tional coefficients for the finite-dimensional resolvent kernel, which improves a 
known result by revealing its rational dependence with respect to the spectral 
variable, as well as with respect to the linear operator. 

In the context of general relativity, the 4-dimensional characteristic formula 
endows a polynomial expression for the volume scalar density. Furthermore, the 
4-dimensional resolvent formula yields a tensor- valued third-degree polynomial 
for the inverse metric, thus avoiding the computacional and formal drawbacks 
of the Neumann series. 



II. Newton's identities 

The fundamental theorem of algebra QXJ| together with Euclid's division al- 
gorithm imply algebraic closureness |2[] of the complex number field; this basic 
result is encompassed without proof by 

Theorem // p(z) = z n + D 1 z n ~ 1 + • • • + D n -\ z + D n is a polynomial with 
complex coefficients then there exist complex numbers Ai, A2, • • • , A n , called 
roots of p, such that p(z) — (z — Xi)(z — A2) ■ ■ ■ (z — X n ). 

As a consequence of the identity principle, the following relations between 
coefficients and roots hold: 

-D x = Ai + A 2 + • • • + A n , 
+D2 = A1A2 + ■ ■ • + A n _iA n , 

(-) k D k = £ A;-products of A's, ^ 

{—) n D n = A1A2 ■ ■ ■ A n . 

The right-hand side of (p]) defines the elementary symmetric functions 
in the variables Ai, A2, . . . , A n . Next to them, the most important symmetric 
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functions are the sums of like powers: 

T x = Ax + A 2 + • • ■ + X n , 
T 2 = A! 2 + A 2 2 + --- + A n 2 , 

k '■■ k k (2) 

Tk = Ai + A 2 + • • • + A n , 

.Dfc and Tjt, besides being symmetric, are homogeneous functions of degree 
k in the variables Ai, A 2 , . . . , A n . We shall derive a set of recursive relations 
connecting D's and T's by which the A's are eliminated from ([I]) and 

To pursuit this goal we need two lemmas. The first one is an improved 
version of the remainder theorem. 

Proposition 1 If p(z) = z n + D\ z n ~ l + ■ ■ ■ + D n -i z + D n is a polynomial 
with complex coefficients Dk and X is a complex number then 

~ P x ^ = zn ~ l + ( A + D i) zn ~ 2 + (A 2 + -Di A + D 2 ) z n ~ 3 + ■■■ 
+ (\ n - 1 + D 1 \ n - 2 + --- + D n _ 1 ). 

Proof; By induction on the degree. For n—1, p(z) — p(X) = (z + D{) — (A + 
Di) = z — A. For generic n, p(z) — p(A) = 

= (z n + D 1 z n ~ l + ■■■ + D n - lZ + D n ) - (A" + D x A"- 1 + • • • + ZViA + D n ) 
= z {z n ~ x + D x z n - 2 + ■■■ + D n ^) - A (A™- 1 + D 1 A"" 2 + • • • + D n _i) 
= z ((z™- 1 + D x z n ~ 2 + ■■■ + D n _ x ) - (A™- 1 + D x A™" 2 + • • • + £>„_i)) + 
+ (2 - A)(A"- 1 + Di X n ' 2 + ■■■ + D n _ x ). 

From £/ie inductive hypothesis the last expression reads 

z(z-X) {z n ~ 2 + (A + Dx) z ra ~ 3 + ■ • • + (A™- 2 + A™~ 3 + ■ ■ ■ + D n _ 2 )) + 
+ (z - AXA"- 1 + D 1 X n - 2 + ■■■ + D n _ x ) 
= (z - X)(z n ' 1 + (A + D x )z n ' 2 + ■■■ + (X n - 2 + DiA n - 3 + ■ ■ ■ + D n _ 2 ) z + 
+ (A"- 1 + J D 1 A"- 2 + --- + J D n _ 1 )). 

The second lemma is a contribution of differential calculus to algebra. 

Proposition 2 If Ai, A 2 , . . . , X n are the complex roots of the nth-degree poly- 
nomial p and p' is its derivative polynomial, then 

p( z ) = L;— j~- 

Proof; From p(z) = d (z — Xi)(z — A 2 ) • • • (z — X n ), it follows 

p'{z) =d {z- X 2 )(z - A 3 ) • • • (z - X n ) +d (z- Ai)(z - A 3 ) • - • (z - A n ) H 

+d0 (,-A0(,-A 2 )...(,-A„_ 1 )^ + J^- + ... + i 44-. 
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Thus, we are able to get 

Theorem 1 (Newton's formulas [|[]) If X x , X 2 , . . . , X n are the complex roots 
of the polynomial p(z) = z n +D x z n ~ 1 +D 2 z n ~ 2 + - ■ -+D n _i z+D n with complex 
coefficients D k and T k = Xi k + X2 k + • • • + X n k then 

T k + D 1 T k _ x + --- + D k _ 1 T 1 + kD k = 0, k = 1, 2, n. 

Proof: Propositions [I] and [| yield 

nz n ~ l + (n- l)D x z n ~ 2 + {n- 2)D 2 z n ~ 3 + ■■■ + D n _ x = 

n ( 

= £ z 71 - 1 + (A fc + B x )z n ' 2 + (X 2 k + D x X k + D 2 )z n ~ 3 + ■■■ 

k=l I 

+(Ar x + D x X n k ~ 2 + ■■■ + D n . 2 X k + L> n -i)} 

= nz 71 - 1 + (Ti + nD x )z n - 2 + (T 2 + D X T X + nD 2 )^~ 3 + • ■ ■ 
+(T n _! + D x T n _ 2 + ■ • • + + nD n _ x ). 

Equating coefficients we obtain 

(n-l)D 1 = T x + nD x , 

(n - 2)D 2 = T 2 + D X T X + nD 2 , 

D n -i = T„-i + D x T n ^ 2 H h D n „ 2 T x + nD n - X , 

from which follow the (n — 1) /irst relations to be shown. Since each X k is a 
root of p it verifies = X r k l + D x X k ~ 1 + ■ • • + D„- X X k + D n . The addition of 
these relations yields 

n 

= E {K + D ^ X V + ■■■ + D n-Ak + D n ) = T n +D x T n _ x +- ■ ■+D n _ 1 T 1 +nD n , 
fe=i 

which is the remaining formula to be proved. 

III. The trace formulas 

The characteristic polynomial || of a complex linear operator T in a finite 
dimensional complex vector space \V is defined by p(z) = det (zl — T), where I 
is the identity operator in \V. The degree of p as a polynomial in the complex 
variable z equals the dimension of \V as a complex vector space; its roots are 
called the characteristic values of T and the set of characteristic values is called 
the spectrum of T. 

These definitions and properties lead to 
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Proposition 3 The trace of a complex linear operator T equals the sum of its 

characteristic values counted with multiplicities. 

Proof: det (zl — T) = (z— \i)(z— A 2 ) • ■ • (z— A n ). By the matrix representation 
T' l j of T we cast the left-hand side as z n -(T\ + T\ + ■ ■ ■ + T n n ) z n ~ l + terms 
of order < n — 2. In the right-hand side we have z n — (Xi + A 2 + ■ ■ ■ + A n ) z n ~ l + 
terms of order < n — 2. Equating coefficients we obtain 

T\ + T\ + • • • + T n n = Ax + A 2 + • • • + A n . 

The roots of unity in the complex field yield the following factorization 
lemma in the associative algebra of complex linear operators. 

Proposition 4 If T and W are commuting complex linear operators and 9 = 
e i2ir/k _ cog 27T _|_ ^ g j n 2n yjjfa a p 0S itwe integer, then 

T k - W k = (T - W)(T - 6 W)(T - 6> 2 W) • • ■ (T - 6' fc " 1 W). 

Proof: As T and W commute we collect terms as 

(T - W)(T - 6 W)(T - # 2 W) • ■ • (T - O^W) = 

= T k + diWT^ 1 + rf 2 W 2 T fc " 2 + ■ • • + 4_iW fe_1 T + d k W k , 

where the complex numbers d\, rf 2 , • • • , dk are the coefficients of the polynomial 
(z-l)(z-6)(z-6 2 )---(z-6 k - 1 ) = z k -l, since (9 j ) k = U^/ k y) k = 1 for 
each j = 0, 1, . . . , k — 1. Therefore d\ = d 2 = . . . = dk-i = and dk = — 1. 

From the multiplicative and homogeneity properties of determinants we get 

Proposition 5 // Ai, A 2 , . . . , A n are the characteristic values of the complex 
linear operator T then Xi k , A 2 fc , . . . , \ n k are the characteristic values of the 
linear operator T k . 

Proof: For the particular case where W = wl, proposition yields 

det (T k - w k l) = det (T - wl) det (T - Owl) ■ ■ ■ det (T - 9 k ^wl) 
= (Ai -«;)••■ (A n - w) (Ax -Bw)--- (A„ - Sw) • • • (Ai - 9 k - 1 w) ■ ■ ■ (A n - 
= (Ai - w){\x -9w)--- (Ai - • • • (A n - u;)(A n - 0io) • • • (A n - 9 k - 1 w) 

= (Ai fc - w k ) ■ ■ ■ (\ n k - w k ). 

Thus det (zl - T k ) = (z- A x fc ) ■■■{z- \ n k ). 

Collecting the former results we are ready to achieve the recursive for- 
mulas enclosed by 
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Theorem 2 If p(z) = z n + D\z n 1 + ■ • • + D n ^\z + D n is the characteristic 
polynomial of the complex linear operator T and = trace (T fc ) then 

T k + Dx r fc _! + • • ■ + D fe _x Ti + = 0, fc = 1, 2, . . . , n. 

Proof: zs a straightforward consequence of theorem [| together with proposi- 
tions [| and [|. 

The key content of theorem shall be stood out by the trace formulas 
statement: the coefficients of the characteristic polynomial of a linear 
operator can be recursively computed as polynomial functions in the 
traces of its successive powers. 

The trace formulas = _Dfc(Xi, T2, . . . , X&) for the coeffi- 

cients of the characteristic polynomial are listed below for k up to 8. 

Di = -T u 

D 2 = 2^1 2 — 2^2' 

d ?j = -\t^ + 1^- \t 3 , 

n 1 71 4 1t2t i Itt 1 It 1 2 It 

4_ 24 Jl — 4 J l i 2+3-'l-t3+8-'2 — 4-M, 

n 1 T 5 1 1 T 3t lfi 271 lT"r2 1 Iff 1 17171 It 

5 — — 120 1 "l" 12 1 2 _ 6 1 3 ~ 8 1 2 + 4 i l J 4+6 2 3_ 5 5 ' 

7~i 1 71 6 1 T 1 1 7 3T 1 1 T 2rp 2 It 2t I71 rp rp 

u %— 720 J 1 — igJl J 2 + Ys Jl J 3 + Ye-tl ~ g i 1 J- 4 ~ g-i 1-i 2-i 3 

1 I7171 1 t3 1 Itt ± 1t2 It 

+ 5 Jli 5~48 i2 +8 i2i4+ 18-'3 -g J 6, 

7~) 1 rp 7 j irp 5t 1 T 4t 1 rp 3rp 2 1 I rp $rp , 1 rp 271 rp (Q\ 

•^7 -- 504oJl + 240 Jl J2 ~ 72 Jl J 3 - 4g J l J 2 + 24 J 1 J 4 + 55 J 1 -* 2 i 3 W 

1 71 2o^ 1 1 71 71 3 I7 rp rp 1 rp rp 2 1 Irp rp 1 rp 2rp 

— iQ-tl -t5 + 48 J l i 2 — g-tl-t2-t4 — ig J l i 3 + g 1 I 1 6 — 24 J 2 1 3 

I 1 rp rp , 1 rp rp 1 

"TTo J 2 J 5 + i2 i 3 i 4 _ 7 J 7) 

n 1 T 8 It 6t i 1 T 5t i 1 T 4t 2 It 4t 

•^8 - 4032oJl - Tiio^l +2 + ^3 + Tg^l ^2 -ggilA 

1 T 3t T i It 3t 1 T 2t 3 i It 2t T l 1 T 2t 2 
-ggJl i 2 i3 + 3pil is-ggJl -*2 +TgJl J 2 i4+3gJl 1 3 

1 ryn 2T I 1 T T 2T 1 rp rp rp 1 rp rp rp , Irp m , 1 7 4 

— T 2 " i l J 6 + 24 1 2 J 3 — Yo J l J 2 J 5 — i2 i l i 3 i 4 + 7 J l-t7 + 384-^2 

1 t 271 l t t 2 j Itt 1 Itt i l t 2 It 

— 32 J 2 J4 - 36^2^3 + if J 2 J 6 + T5 -* 3-t 5 + 32 1 4 ~ s 1 ^ 

The characteristic formulas for dimensions up to four are written down: 

det (zl — T) = z-T u 

det(zI-T) = z 2 -T l z+ l -{T 1 2 -T 2 ), 

det(zI-T) = Z ^-T 1 z 2 + -\(T 1 2 -T 2 )z--\(T 1 3 -3T 1 T 2 + 2T 3 ), (4) 

2 o 

det {zl - T) = z 4 - Tiz 3 + i(T x 2 - T 2 )z 2 - i(Tx 3 - 3TxT 2 + 2T 3 )^ + 

2 o 

+ 777( r i 4 - 671 2 T 2 + 87^ + 3T 2 2 - 6T 4 ). 
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IV. Null identities 



The polynomial p(z) = d§z n + d\z n 1 + ■ • • + d n -\Z + d n with complex 
coefficients dk is said to annihilate the complex linear operator T if p(T) = 
d T n + diT n_1 H h d n _iT + d n I = O, the identically null operator. 

The knowledge of the principal ideal || of polynomials that annihilate a 
linear operator T is essential to attain computational skill in the associative 
algebra generated by T. To ensure this goal we state without proof one of the 
fundamental results in linear algebra. 

Theorem 3 (Cayley-Hamilton [[7]]) The characteristic polynomial of a lin- 
ear operator annihilates it. 

Joining Cayley-Hamilton's theorem with trace formulas statement we con- 
clude that for each finite dimension of the underlying vector space 
there is a fundamental null identity in the associative algebra of lin- 
ear operators. For instance, in dimensions up to four, the characteristic 
formulas (^) yield the following null identities: 

T — Til = O, 

T 2 — 7\T + 7j(7\ 2 - T 2 )I = O, 

T 3_ TiT 2 + 1 ^2 _ _1 ^3 _ 3TiT2 + 2Ta)I = q 

T 4 -TiT 3 +^(Ti 2 - T 2 )T 2 -^(T! 3 - 3T ± T 2 + 2T 3 )T+ 

+ ^(T X 4 - 6Ti 2 T 2 + 8TxT 3 + 3T 2 2 - 6T 4 )I = O. 

V. The finite-dimensional resolvent kernel 

The resolvent of a complex linear operator T is the operator- valued function 
R defined by R(z) = (zl - T) . It is a well-known fact that R is an operator- 
valued analytic function outside the spectrum of T. We shall refine such a 
result showing that in the finite dimensional case R is an operator-valued 
rational function completely tied by a finite set of complex-valued rational 
functions, whose coefficients are exactly the trace formulas. 

In this vein we need an improved version of the remainder theorem in the 
associative algebra of complex linear operators, whose content is the extension 
of proposition [l] to operator- valued polynomials. 

Proposition 6 If p(w) = w n + Di%v n ~ l + • • • + D n _\w + D n is a polynomial 
with complex coefficients and T and W are commuting complex linear 
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operators then 
p(W) -p(T) = (W-T) 



T"- 1 + (W + DJ)T 1 - 2 + (W 2 + DxW + D 2 I)T r ' 
+ ... + (W"" 1 + DxW"^ 2 + • • • + D n _iI)I . 



Proof: By induction on the degree. For n = 1, p(W) — p(T) = (W + Dil) — 
(T + Dil) = W - T. For generic n, p(W) - p(T) = 



W(W n " i + D^ 1 - 2 + 
-T(T™" 1 + F^T™" 2 + 



(W-T)(W n_1 + £>iW 



m-2 



D n _xl) + F>J+ 
D n _xl) - D n I 
!-••• + A»-iI)+ 



+T 



(rr\n— 1 



n-2 



DlT n-2 + 



+ • • • + FW) + 



From £/ie inductive hypothesis the last expression reads 

(w - t) (w n_1 + r^w™- 2 + • ■ ■ + £>„_il) 



+T(W - T) 



(W-T) 



T™- 2 + (W + F^T™" 3 + ■ • • 
+ (W"- 2 + F^W^ 3 + ■■■ + F> n _ 2 I)I 
T"" 1 + (W + Dxl)^ 2 + • 
+ ■ • • + F> n _ 2 I)T + (W"- 1 



• + (W n ~ 2 + DxW n - 3 + 
F>iW n ~ 2 + ■ • • + F> n _iI)I 



We are ready to achieve a rational formula for the resolvent kernel. 



Theorem 4 If p(w) 



w 



Diw n + • • • + D n _iw + D n is the characteristic 



polynomial of the complex linear operator T and z is any complex number not 
belonging to the spectrum of T then 



[zl 



m- 1 



p(z) p(z) 
. z n ~ l + D x z n - 2 



Z + Djrpn-2 , Z 2 + D 1 Z + D 2 r r n-3 



+ 



pjz) 



+ 



p{z) 



I. 



Proof: It is enough to consider Cayley-Hamilton's theorem for the linear op- 
erator T and set W = zl in proposition [6|/ notice that in such case W fe + 
DjW^ 1 + ■■■ + D k l = (z k + D 1 z k ~ 1 + ■■■ + D k )I for each k = 1, 2, . . . , n. 

The joint content of theorem [| and the trace formulas statement shall be 
pointed out as : for each finite dimension of the underlying vector space 
there is a fundamental rational formula for the resolvent kernel of a 
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linear operator. For instance, in dimensions up to four, the trace formulas 
(0) provide the following resolvent formulas: 



zl- 


T) 


-^(z-T,)- 1 !, 






zl- 


T) 


- 1 =(z 2 -T l z + 






zl- 


T) 


1 = (** —Tiz 2 + 




-T 2 )z- 



(zl - TV 1 



t + (z — r x )i 

o 



-1 



T 2 + 



+(z - Ti)T + [z 2 - T lZ + -(71 2 - T 2 )J I 

- T lZ 3 + U^ 2 - T 2 )z 2 - i(Tx 3 - 3T X T 2 + 2T 3 )z+ 
2 o 



+— (Tx 4 - 6T! 2 T 2 + 8T!T 3 + 3T 2 2 - 6T A ] 



-i 



T 3 + 



+ (z - Tx)T 2 + [z 2 - T lZ + -(71 2 - T 2 )J T + (z 3 + 
-T^ 2 + i(T x 2 - T 2 )z - 1 -{T^ - 2T X T 2 + 2T 3 )) I 



(5) 



VI. Applications to general relativity 

In general relativity the fundamental physical object is an effective geom- 
etry, mathematically represented by a second-rank covariant tensor field g, 
defined throughout a suitable 4-dimensional manifold. In a local coordinate 
system x = (x a ) the metric tensor can be written as g = g^dx^ ® dx u . 

Investigations on general relativity are frequently carried out under the 

o o 

assumption that there exists some background geometry g= 9 ^dx^ <8> dx v . 
The metric properties to be assumed on g vary depending on the gravitational 
scenario. In order to perform calculations it suffices to set g a Ricci-flat metric. 
However, to interpret the results as physically meaningful it is generally agreed 
that one should require g to be a flat metric. Even this case is sometimes 
thought to be too broad, as some authors claim to set harmonic coordinates 
[U or even cartesian coordinates 0. We do not take into account here any 
suplementary conditions on the background geometry g. 

The effective and background geometries are related by a tensor field h = 
h^dx 11 ® dx u by a connecting equation, the generally accepted form of which 
being g = g + h. The scalar density y/—gdx 4 associated with g requires us to 
compute det (g) = det (g +h). This can be achieved by means of 
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Proposition 7 If g =g + h, H =g 1 h, H k = trace (H fc ) and i/ie underlying 
manifold is 4- dimensional then 



det 



det 



1 + H x +\{H 1 2 - H 2 ) + ^(H^ - 3#i# 2 + 2# 3 ) + 

+ ^(i/x 4 - %H X 2 H 2 + + 3# 2 2 - 6# 4 ). 



Proof: From g =g +h =g (I + H) it follows that det (g) = det (g) det (I + H). 

Now it is enough to consider the 4-dimensional characteristic formula in ffl) 
with z = 1 and T = -H; notice that T k = (-) k H k implies T k = (-) k H k . 

The Levi-Civita connection associated with g requires us to compute g _1 = 
(g +hj = (I + H) 1 g _1 . The known explicit form is given by the Neumann 



series 



10 



(I + H)" 1 = I- H + H 2 -H 3 + 
In local coordinates it reads 

<T = (V/3- 9*"* haf>+ °9^h aX °g Xe h e/3 - °g^h aX °g^K p g^Kp + - ■ • ^ " 11 

o o „ o 

where g^ u and 9^ are well-defined by g m g a u = ^ u =g^ p 9 p u . 

Besides convergence requirements on the above series, we stress that such 
expression leads to technical difficulties when developing the Lagrangian varia- 
tional formalism. These problems were already dealt with in the literature, the 
proposed solution being to modify the above form of the connecting equation 
[PH . We show how to overcome such drawbacks by means of 

Proposition 8 If g =g + h, H =g~ 1 h, H k = trace (H fc ) and the underlying 
manifold is 4-dimensional then 

g- 1 = (l + H 1 + ^{H, 2 - H 2 ) + \W - SH X H 2 + 2H 3 ) + 

+^(#i 4 - m x 2 H 2 + 8#i# 3 + 3H 2 2 - QH A )y l - H 3 + 

+ (1 + #i)H 2 - (l + H x + \{H\ 2 - H 2 )) H + (l + 

+#! + ^(H, 2 - H 2 ) +^(#i 3 - 3i7xF 2 + 2H 3 j) I 



" -1 
g • 



Proof; From g =g + h =g (I + H) it follows that g 1 = (I + H) 1 g 1 . Now 

it is enough to consider the 4-dimensional resolvent formula in (|3p with z = 1 
and T = -H; notice that T k = (-) fc H fc implies T k = (-) k H k . 
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